Ultracold neutral plasma expansion in two dimensions 
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We extend an isothermal thermal model of ultracold neutral plasma expansion to systems with- 
out spherical symmetry, and use this model to interpret new fluorescence measurements on these 
plasmas. By assuming a self-similar expansion, it is possible to solve the fluid equations analytically 
and to include velocity effects to predict the fluorescence signals. In spite of the simplicity of this 
approach, the model reproduces the major features of the experimental data. 
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Ultracold plasmas are produced from photo-ionized 
laser-cooled gases pi 0, Q • In these laboratory plasmas, 
it is possible to study the kinetics and thermodynamics 
of multi-component, strongly-interacting Coulomb sys- 
tems. These systems are characterized by the ratio of 
the nearest-neighbor Coulomb energy to the average ki- 
netic energy, denoted as T = (e 2 /ineod)/ (kbT), with d 
being the interparticle spacing. 

A new class of ultracold plasma experiments has re- 
cently become available in which it is possible to spec- 
troscopically study the plasma ions 0, 0, @ • These plas- 
mas are made using alkaline-earth atoms, because the 
resonance transition wavelengths of the ions are readily 
generated using standard laser methods. The spatially- 
resolved time evolution of the plasma ion temperature 
and density can be measured using absorption and fluo- 
rescence techniques. 

These ultracold neutral plasmas are not trapped, al- 
though efforts are underway in a few laboratories to trap 
them. The untrapped plasmas freely expand, and as they 
expand the density and temperature change radically. 
Processes of recombination, collisional and thermal ion- 
ization, radiative cascade, and adiabatic and evaporative 
cooling all play important roles in how the system evolves 
and equilibrates. 

A variety of models have been used to investigate the 
properties of these plasmas H H El El 13 El El 
One particularly simple isothermal fluid model |8( has 
been surprisingly successful in predicting the general fea- 
tures of these plasmas 0, H 0. In this paper we ex- 
tend this model from the spherically-symmetric Gaus- 
sian plasma distributions to Gaussian distributions with 
elliptical symmetry. 

The elliptical symmetry has important experimental 
advantages. In such systems the plasma expands primar- 
ily in two dimensions. The practical advantage is that the 
density falls more slowly than in the three dimensional 
case, making it possible to study the plasma for longer 
times. The Doppler-shift due to the directed expansion 



of the plasma is also suppressed. It should therefore 
be possible to study plasma oscillations and heating ef- 
fects for greater time periods before these oscillations are 
masked by the directed expansion of the plasma. Finally, 
if the plasmas are generated from a density-limited neu- 
tral atom trap, the elongated symmetry allows a greater 
number of atoms to be trapped initially, corresponding 
to a greater column density of plasma ions. This directly 
increases the visibility of fluorescence and absorption sig- 
nals. 



I. ISOTHERMAL FLUID MODEL 

An isothermal fluid model has been presented in the 
literature 0, • It successfully reproduces most of the 
major features of recent experimental work. This model 
was motivated by trends observed in more sophisticated 
treatments. The basic ideas of the model will be reviewed 
here, and an extension to the case of a Gaussian distri- 
bution with elliptical symmetry will be presented. 

The initial ion density distribution is proportional to 
the Gaussian distribution of the neutral atom cloud from 
which the plasma is created. For a spherically sym- 
metric cloud, the initial distribution can be written as 
n(r) = uq exp(— (3r 2 ). Because the electrons thermalize 
much faster than the ions, in this model we take the ini- 
tial electron density distribution to be the thermal equi- 
librium distribution, given by the Boltzmann factor: 



n e (r) = n 0e exp 



eV(r) 
k P ,T 



(1) 



The lowest temperature plasmas are nearly charge- 
neutral, and the electron density distribution is approxi- 
mately equal to the ion density. In this limit, it is shown 
m HI that for a spherically symmetric plasma that to 
within an arbitrary additive constant the electrical po- 
tential energy can be written as 
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eV(r) = k B Tbx[n{r)/n 0e ] = -k B Tf3r 2 . (2) 

The force is the negative gradient of this potential en- 
ergy. It is manifestly radial, spherically symmetric, and 
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linearly proportional to the radial coordinate r, measured 
from the center of the plasma. The velocity, which is the 
time integral of the acceleration, is also linearly propor- 
tional to r. The consequence is that if the distribution is 
Gaussian initially, it will remain Gaussian at all times in 
the expansion. 

For the case of non-spherical symmetry, the approach 
is more or less the same, although the isothermal nature 
of the plasma has a more restricted meaning. We will 
take the initial ion density distribution to be Gaussian, 
symmetric in the x — y plane, and initially elongated in 
the z direction: 



Equation[S]can be re-written (dropping the subscripts) 



n(r, i) 



_3/2 



exp Hz 2 + y 2 )/3i(t)-z 2 /3 2 (t)] , 



(3) 

The initial conditions are /?i(0) = <j 2 and /^(O) = 
a 2 aQ 2 , and the parameter a defines the elipticity of the 
system. The plasma fluid equations for our system are 
written 



+ V • (nv) = 



dn 

dv _ 

— + (v ■ \7)v = a. 



(4) 
(5) 



Following the derivation used in the case of spherical 
symmetry, the velocity is taken to be v = 'yx(t)(xx+yy) + 
j2(t)zz. Inserting this and the density profile of Eq. 03 
into Eq. Ogives 



2(3 2 (x 2 p 1 + y 2 /?! - 1) (ft + 2/3 l7l ) • • • 

+ 4/?i (z 2 /3 2 - 1) (p 2 + 2/3 272 ) = 0. (6) 

Because x, y, and z are independent variables, the only 
non-trivial solution of this equation is 



7 = 



(7) 



where we have dropped the subscripts because all com- 
ponents have this same form of solution. 

Solving Eq. [5]requires a little more care. It is straight- 
forward to write down the acceleration vector follow- 
ing the derivation of the spherically symmetric case. 
However, the solution requires that the temperature be 
isothermal in a given dimension, but anisotropic in space. 
This condition allows the density distribution to reduce 
to the proper form in the limiting case of a plasma in- 
finitely long in the z dimension. This decouling requires 
energy to be conserved separately in the x — y plane and 
in the z dimension. The plasma equations for these two 
spaces are now exactly identical and completely separa- 
ble. 



as 



7 + 7 2 = 2k b T(t)P(t)/m, 
and the conservation of energy is 



T(0) = T(t) 



m 7 
2k~ b J 



(8) 



(9) 



where we have neglected the energy due to electron-ion 
recombination. Equations 00 and are exactly iden- 
tical to Eq. 2 of Ref. ||. Using Eqs. and 03 we now 
have both T and 7 in terms of (3. Inserting this into Eq. 
IH1 gives 



0i u > 



(10) 



where we have made the substitution v 2 = 2ki,T(0)/m. 
The solution to this equation is /3f 1 = a\ + ct + v 2 t 2 , 
where c is an integration constant. The constant c must 
be equal to zero to meet the condition that the ion ve- 
locity is initially zero at t = pfl. The time evolution of 
the density and velocity functions can now be written as 



n(r,t) = 



v(f, t) 



iVexp [-(x 2 + y 2 )/a 2 {t) (z 2 )/a 2 (t)} 

^/ 2 a 2 {t)a 2 {t) 1 > 



vlt 



a((t) a 2 {ty a 2 (t) 



a 2 (t) 
4(t) 



1 1 2,2 
a +v e t 

a 2 ol+v 2 A 2 . 



(12) 

(13) 
(14) 



We note that for two-dimensional planar Gaussian 
charge distributions, a closed- form expression for the 
electric field has been derived 0, 0|. If such a com- 
pact analytical solution could be written in the three- 
dimensional case, it would remove the decoupling con- 
straint that we imposed in order to solve the plasma 
equations. However, such a solution is not readily ap- 
parent. 



II. FLUORESCENCE SIGNAL MODEL 

The geometry of our measurements is shown in Fig. 
In the experiment, the probe laser beam is spatially fil- 
tered and focused into the plasma with a confocal beam 
parameter that is long compared to all plasma dimen- 
sions. The position offset of this probe laser relative to 
the plasma is denoted by the parameter a. After the 
plasma is created, the number of atoms in the column 
defined by the probe laser beam changes dramatically. 
When the laser beam passes through the center of the 
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FIG. 1: A schematic diagram of the fluorescence measure- 
ments. Left — The plasma is initially at a high density, and 
elongated in the ^-dimension (labeled region 1). The probe 
laser is represented by the solid dark arrow parallel to the z 
axis, displaced by a variable distance a. At later times the 
plasma expands slowly in z, but quickly in the x — y plane. 
Regions 2 and 3 represent the plasma size at later and later 
times. As the plasma expands, the number of atoms in the 
column defined by the probe laser changes. Right — a partial 
level diagram of Ca + . The probe laser is locked to the 397 nm 
resonance transition, and we measure 397 nm light scattered 
by plasma ions. 



plasma (a — 0), the number of ions in the beam mono- 
tonically decreases. However, if the laser is outside of 
the initial plasma distribution, as shown in Fig. ^ the 
number first increases as the plasma moves into the laser 
beam, and then decreases as the plasma disperses. 

The fluorescence signal depends on both the number 
of ions in the column defined by the probe laser beam 
(Gaussian beam profile with 1/e 2 radius w) and the ve- 
locity distribution of the ions in the plasma. Because the 
laser has a narrow bandwidth, atoms moving at veloci- 
ties greater than ~9 m/s are Doppler-shifted out of reso- 
nance. In this section we will use the results of the previ- 
ous section to derive an expression for how the plasma ion 
fluorescence signal should change with time for different 
values of the offset parameter a. 

The fluorescence signal s(t) is proportional to the ab- 
sorption of the probe laser beam. Using Beer's law and 
the standard approximation of small optical depth, s(t) 
can be written as 



s(t) oc / n{x,y,z,t)f\(x - a,y)a{u - v ), (15) 
./Vol 

where /i is the spatial profile of the probe laser beam, and 
a is the absorption cross section as a function of v — 
the difference between the laser frequency and the atomic 
resonance frequency. Removing the limitation of small 
optical depth is trivial. We can simplify this expression 
by setting the laser frequency equal to Uq, and recogniz- 
ing v — vq — v/X as the Doppler shift due to the velocity 



v of the atoms, where A is the optical wavelength of the 
transition. Ea uat ion 1121 gives the position-dependent ve- 
locity of the ions. We define a length £ = T\a 2 /2v 2 t, 
where T = 1/2ttt is natural width of the transition, r is 
the lifetime of the transitions's upper state, and use it to 
write the absorption profiles, 



ct/ctq = 



(z/t) 2 



exp 



-2(v th /2TX) 2 (z/£) 2 



Lorentzian, 



Gaussian, 



(16) 

where Co is the absorption cross-section on resonance and 
Vth is the rms velocity of a thermal distribution. The true 
absorption lineshape is better represented by a Voigt pro- 
file. However, as the Voigt profile can be approximated 
by a linear combination of the Lorentzian and Gaussian 
line profiles |17| , we will only write down these two limit- 
ing forms. Power broadening of the line can be included 
in a straightforward manner |l8j | . 

We take the plasma density profile from Eq. ^] and 
write the spatial profile of the probe laser beam as 



h = exp [-2(a; - a) 2 /w 2 - 2y 2 /w 2 ] , (17) 

which corresponds to the geometry represented in Fig. ^ 
Performing the integration in Eq. 1151 gives. 



rj erfc(?7) exp(?7 2 ) 



Lorentzian, 



s(t) oc — exp(-£) 

1 + {2TX/v th i]) 2 Gaussian, 

(18) 

where £ = 2a 2 /(w 2 + 2a 2 ) and r\ = lju2- This expression 
for the Lorentzian lineshape is proportional to Eq. 5 of 
Ref. H. 



III. ULTRACOLD CALCIUM PLASMAS 

We create ultracold neutral plasmas by photoioniz- 
ing laser-cooled calcium atoms in a magneto-optical trap 
(MOT). The calcium MOT is formed in the usual way 
by three pairs of counter-propagating laser beams that 
intersect at right angles in the center of a magnetic 
quadrupole field 0]. The 423 nm laser light required 
for the calcium MOT is generated by frequency-doubling 
an infrared laser in periodically-poled KTP (PPKTP), 
and has been described previously [2(j. A diode laser 
master-oscillator-power-amplifier (MOPA) system deliv- 
ers 300 mW single frequency at 846 nm, as shown in Fig. 
121 This laser is phase-locked to a build-up cavity using 
the Pound-Drever-Hall technique [2lJ , giving a power en- 
hancement of 30. A 20mm long PPKTP crystal in the 
small waist of the build-up cavity is used to generate 
typically 45 mW output power at 423 nm [2^, [2jj . 
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FIG. 2: A schematic drawing of the MOT laser system and 
frequency stabilization electronics used in these experiments. 



The laser is further stabilized by locking the 423 nm 
light to the calcium resonance transition using saturated 
absorption spectroscopy in a calcium vapor cell [24| . Our 
vapor cell differs from Ref. [24| in that it has a stainless 
steel body with conflat metal seals and windows and a 
valve. An acousto-optic modulator (AOM) in one arm 
of the saturated absorption laser beams shifts the laser 
frequency so that the laser beam sent to the MOT is 
35 MHz (one natural linewidth) below the atomic res- 
onance. We also use the AOM to chop this beam and 
use a lock-in amplifier to eliminate the Doppler back- 
ground in the saturated absorption signal. Because the 
846 nm laser is already locked to the frequency-doubling 
cavity, the feedback from this second lock circuit servos 
the frequency-doubling cavity length. 

The trap is loaded from a thermal beam of calcium 
atoms that passes through the center of the MOT. The 
thermal beam is formed by heating calcium in a stainless 
steel oven to 650° C. The beam is weakly collimated by 
a 1mm diameter, 10mm long aperture in the oven wall. 
The distance between the oven and the MOT is approx- 
imately 10 cm. As the beam passes through the MOT, 
the slowest atoms in the velocity distribution are cooled 
and trapped. An additional red-detuned (140 MHz, or 
four times the natural linewidth) laser beam counter- 
propagates the calcium atomic beam, significantly en- 
hancing the MOT's capture efficiency. To prevent op- 
tical pumping into metastable dark states we also em- 
ploy a diode laser at 672 nm. The density profile of the 
MOT has an asymmetric Gaussian profile and is well- 
represented by Eq. 1111 with the peak density equal to 
4 x 10 9 cm" 3 , t7o = 0.5 mm, and a = 2.5. 

We photo-ionize the atoms in the MOT using a two- 
color, two-photon ionization process. A portion of the 
846 nm diode laser radiation from the MOT laser is 
pulse-amplified in a pair of YAG-pumped dye cells and 
frequency doubled. This produces a 3 ns-duration laser 
pulse at 423 nm with a pulse energy around 1 /zJ. This 
laser pulse passes through the MOT, and its peak in- 



FIG. 3: A schematic diagram of the DAVLL cell. Upper — A 
drawing of the hollow. cathode discharge cell, approximately 
to scale. The cathode is 10 cm long, with a 1 cm diameter 
bore. The operating pressure is between 100 and 200 mTorr 
of krypton. Lower — Electrical connections to the discharge 
cell. For increased current stability, a 2kf2 ballast resistor is 
connected in series with the discharge. 



tensity is a few thousand times greater than the satura- 
tion intensity. A second YAG-pumped dye laser at 390 
nm counter-propagates the 423 nm pulse and excites the 
MOT atoms to low-energy states in the region of the ion- 
ization potential. We photo- ionize 85-90% of the ground- 
state atoms in the MOT. The minimum initial electron 
temperature is limited by the bandwidth of the 390 nm 
laser to about 1 K. 

Ions in the plasma scatter light from a probe laser 
beam tuned to the Ca II 2 <Si/2 — 2 -P1/2 transition at 397 
nm. The probe laser is generated by a grating-stabilized 
violet diode laser. This laser, as well as the 672 nm laser 
used in the neutral atom trap, is locked to the calcium 
ion transition using the DAVLL technique t 25| in a large- 
bore, low-pressure hollow cathode discharge of our own 
design (see Fig. |3J|. 

The probe laser is spatially filtered. The typical probe 
laser intensity is a few hundred fiW focused to a Gaus- 
sian waist of 130 /im in the MOT. We average repeated 
measurements of the scattered 397 nm radiation with the 
probe laser in a given position, denoted by the param- 
eter a in Fig. ^ This produces a time-resolved signal 
proportional to the number of atoms resonant with the 
probe beam in a particular column of the plasma. By 
translating a mirror just outside the MOT chamber, we 
scan the probe laser across the ion cloud. In this manner 
we obtain temporal and spatial information about the 
plasma expansion. 



IV. COMPARING THE MODEL TO THE DATA 

One comparison of the isothermal model with exper- 
imental data was presented in Ref. m. In that work, 
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the initial electron energy of the plasma, and therefore 
the expansion velocity v e , was fixed, and the parameter 
a was varied from to 4uo. 

In the following, we present a complementary compar- 
ison. The solid line in Fig. 21 shows ion fluorescence 
signal with the probe laser tuned to the ion resonance 
frequency, and with the probe laser beam propagating 
along the z axis (a = in Fig. ^| for a range of initial 
electron energies. 

As discussed in Ref. 0, the early rise in the fluores- 
cence signal shows the increasing number of plasma ions 
in the 2 P\/2 level. This feature is easily explained in 
terms of the classic Rabi two-level atom with damping. 
The probe laser intensity is in the range of 5 to 10 times 
the resonance Rabi frequency. Our numerical integration 
of the optical Bloch equations mimics the approximately 
20 ns rise time observed in the experimental data, and 
shows that the Gaussian spatial profile of the probe laser 
beam washes out subsequent oscillations in the excited 
state fraction. Following the initial rise, the fluorescence 
decays. At approximately t = 50 ns the decay slows down 
due to correlation-induced heating in the plasma 0, . 

These two processes are not included in the isother- 
mal model. We therefore begin the comparison of the 
model with the data at time t = 10~ 7 s. The model 
is normalized to the data and fit using v e as the only 
fit parameter in an un-weighted least-squares fitting rou- 
tine. The model uses the Lorentzian lineshape in Eq. 
IT51 The justification for using this lineshape arises from 
the analysis presented in Ref. 0. In Fig. 2 of that 
reference, the t = velocity is 6 m/s. This is the rms 
velocity of a Boltzmann distribution, and it is due al- 
most entirely to correlation-induced heating. This veloc- 
ity width gives a Doppler width smaller than the natu- 
ral line width, and the Voigt profile is close to a pure 
Lorentzian. Furthermore, as the plasma evolves, the ion 
temperature falls due to the adiabatic expansion. The 
time scale for this is ~ (Jo/v e = 2 /is. This is shorter 
than the time scale for heating the ions by collisions with 
the electrons, ~ mQ a /m e ujp = 50 fis. It is therefore not 
surprising that the Gaussian lineshape gives a poorer fit 
to the data. In the fitting procedure, the Gaussian line- 
shape produces expansion velocities that do not grow as 
the initial electron energy increases. The figure shows 
that the Lorentzian model describes the data well over a 
few orders of magnitude in time. The velocities extracted 
from these data are shown in Fig. |SJ 

The differences between the model and the signal are 
not negligible. For all initial electron energies, the model 
is slightly too low at t — 300 ns, and too high at t > 1 fis. 
The data are not corrected for optical pumping into the 
D— states, which we measure to have a time constant of 
about 5 lis. The differences between the model and the 
fluorescence signal could indicate internal heating pro- 
cesses that manifest themselves in the ion velocity on 
the few hundred ns time scale. They could also indicate 
ions that appear in the plasma from Rydberg states at 
late times or collective plasma density variations. These 
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FIG. 4: Relative ion fluorescence signal. The probe laser is 
tuned to resonance and propagates through the center of the 
plasma (a — 0). The electron energy, E e , is measured rela- 
tive to the ionization limit. The solid line is the fluorescence 
signal, the dotted line is a fit using the Lorentzian lineshape 
in Eq. 1181 The model is normalized to the fluorescence signal 
at t = 10 -7 s, and the single fit parameter is the expansion 
velocity, v e . In the top panel, the plasma is spontaneously 
generated after exciting the atoms to a Rydberg state with 
n ~ 35.7. 



differences can be studied in future work. 

We also compare the fluorescence signal and the model 
over a range of initial electron energies with the probe 
laser beam shifted to a = 1.1mm = 2.2oo- In this ar- 
rangement, the ion fluorescence signal is initially small, 
and grows as ions move into the probe laser beam. Typi- 
cal data are shown in Fig. OD For these data, the model is 
fit to the t > 6 /is signal using a least-squares procedure, 
with v e as the fit parameter. The velocity extracted from 
this fit is shown in Fig. It is possible to fit the data 
using v e so that the peak in the model coincides with the 
peak of the fluorescence signal. These velocities are also 
plotted in the figure. 

For these data, the signal in the first /xs is small com- 
pared to the model. Moreover, there are variations in 
the fluorescence signal that do not appear in the model. 
This suggests that at the edges of the plasma expansion, 
the density distribution is distinctly non-Gaussian, even 
at early times before any ion motion is possible. It also 
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time (s) 

FIG. 5: Relative ion fluorescence signal. The probe laser is 
tuned to resonance and propagates at a distance a — 1.1mm = 
2.2<7o relative to the center of the initial distribution. The 
solid line is the data, and the dotted line is the fit using the 
Lorentzian lineshape in Eq. 1181 The model is normalized to 
the peak of the data. 

appears, as pointed out in Ref. [6j, that the Gaussian 
density profile is recovered at late times for all initial 
electron energies. 

V. CONCLUSION 

We present an extension of the isothermal plasma ex- 
pansion model of Refs. 0, 0] for quasi-two-dimensional 
geometry. We include velocity effects and predict a flu- 
orescence or absorption signal vs. time for given initial 
conditions. The model matches the correct order of mag- 
nitude and general features of the fluorescence signal. 
Some discrepancies are pointed out, which can be studied 
in future work. 

Making the plasmas more ideally two-dimensional 
should improve applicability of the model, and further 
suppress effects due to expansion in the long dimension. 
Increasing the range of initial electron temperatures and 
plasma densities can probe interesting regions of phase 



space where differences between the data and the model 
are likely to be more pronounced. 




FIG. 6: Expansion velocity extracted from Figs. 0] and [£] 
The triangles (a) represent the velocity extracted from Fig. 
^] and the circles (•) from Fig. The stars represent 
the velocity extracted from centering the peak of the model 
on the peak of the data in Fig. The line is a fit of the data 
with E e > 150 cm -1 . The fitted result is E e — 3.0mv^, where 
m is the mass of the calcium ion. 



As an example, it should be possible to extend this 
model to include effects due to electron-ion recombina- 
tion at early times. For strongly-coupled neutral plas- 
mas, the three-body recombination rate should be on the 
order of the plasma frequency. Using the high-sensitivity 
and fast time-response of fluorescence spectroscopy, it 
should be possible to directly measure spectroscopic re- 
combination signatures in low-temperature, low-density 
plamsas, where the predicted three-body recombination 
rate is greater than the plasma frequency. For exam- 
ple, reducing the plasma density to 10 6 cm -3 will make 
the recombination time ~ 1/lo p = 20 ns. Increasing the 
probe laser intensity to 100 times the saturation intensity 
will shorten the early rise time of the fluorescence signal 
to around 10 ns. Optical pumping time will be compara- 
ble to the correlation-induced heating time, around 200 
ns. 
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